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Abstract. A quasi-Lie scheme is a geometric structure that provides t-dependent changes of 
variables transforming members of an associated family of systems of first-order differential 
equations into members of the same family. In this note we introduce two quasi-Lie schemes 
for studying second-order Gambier equations in a geometric way. This allows us to study 
the transformation of these equations into simpler canonical forms, which solves a gap in the 
previous literature, and other relevant differential equations, which leads to derive new con- 
stants of motion for families of second-order Gambier equations. Additionally, we describe 
general solutions of certain second-order Gambier equations in terms of particular solutions 
of Riccati equations, linear systems, and t-dependent frequency harmonic oscillators. 
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1 Introduction 

Apart from their inherent mathematical interest, differential equations are important due to 
their use in all branches of science [41, 50]. This strongly motivates their analysis as a means 
to study the problems they model. A remarkable approach to differential equations is given 
by geometric methods [47], which have resulted in powerful techniques such as Lax pairs, Lie 
symmetries, and others [48, 49]. 

A particular class of systems of ordinary differential equations that have been drawing some 
attention in recent years are the so-called Lie systems [1, 9, 33, 46, 54, 56]. Lie systems form 
a class of systems of first-order differential equations possessing a superposition rule, i.e. a func- 
tion that enables us to write the general solution of a first-order system of differential equations 
in terms of a generic collection of particular solutions and some constants to be related to initial 
conditions [7, 18]. 

The theory of Lie systems furnishes many geometric methods for studying these systems [6, 
17, 19, 21, 22, 26, 27, 56]. For instance, superposition rules can be employed to simplify the use 
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of numerical techniques for solving differential equations [56] , and the theory of reduction of Lie 
systems reduces the integration of Lie systems on Lie groups to solving Lie systems on simple 
Lie groups [6, 19]. 

The classification of systems admitting a superposition rule is due to Lie. His result, the 
nowadays called Lie-Scheffers theorem, states that a system admits a superposition rule if and 
only if it describes the integral curves of a t-dependent vector field taking values in a finite- 
dimensional Lie algebra of vector fields. The existence of such Lie algebras on M and M 2 was 
analysed by Lie in his famous work [45]. More recently, the topic was revisited by Olver and 
coworkers [28] , who clarified a number of details that were not properly described in the previous 
literature. 

Despite their interesting properties, Lie systems have a relevant drawback: there exist just 
a few Lie systems of broad interest [36] . Indeed, the Lie-Scheffers theorem and, more specifically, 
the classification of finite-dimensional Lie algebras of vector fields on low dimensional mani- 
folds [28, 45] clearly show that that being a Lie system is the exception rather than the rule. 
This has led to generalise the theory of Lie systems so as to tackle a larger family of remark- 
able systems [2, 3, 9, 36]. In particular, we henceforth focus on the so-called quasi-Lie schemes. 
These recently devised structures [15, 20] have been found quite successful in investigating trans- 
formation and integrability properties of differential equations, e.g. Abel equations, dissipative 
Milne-Pinney equations, second-order Riccati equations, and others [9]. In addition, the ob- 
tained results are useful so as to research on the physical and mathematical problems described 
through these equations. 

In this work, we study the second-order Gambier equations by means of the theory of quasi- 
Lie schemes. We provide two new quasi-Lie schemes. Their associated groups [15] give rise to 
groups of t-dependent changes of variables, which are used to transform second-order Gambier 
equations into another ones. Such groups allow us to explain in a geometric way the existence 
of certain transformations reducing a quite general subclass of second-order Gambier equations 
into simpler ones. Our approach provides a better understanding of a result pointed out in [30]. 
As a byproduct, we show that the procedure given in the latter work does not apply to every 
second-order Gambier equation, which solves a gap performed in there. 

We provide conditions for second-order Gambier equations, written as first-order systems, 
to be mapped into Lie systems via t-dependent changes of variables induced by our quasi-Lie 
schemes. This is employed to determine families of Gambier equations which can be transformed 
into second-order Riccati equations [31], Kummer-Schwarz equations [4, 5] and Milne-Pinney 
equations [25] . These results are employed to derive, as far as we know, new constants of motion 
for certain second-order Gambier equations. Moreover, the description of their general solutions 
in terms of particular solutions of t-dependent frequency harmonic oscillators, linear systems, 
or Riccati equations is provided [9, 29]. 

The structure of our paper goes as follows. Section 2 addresses the description of the fun- 
damental notions to be employed throughout our work. Section 3 describes a new quasi-Lie 
scheme for studying second-order Gambier equations. In Section 4 this quasi-Lie scheme is used 
to analyse the reduction of second-order Gambier equations to a simpler canonical form [30]. 
By using the theory of quasi-Lie systems, we determine in Section 5 a family of second-order 
Gambier equations that can be mapped into second-order Kummer-Schwarz equations. The 
investigation of constants of motion for some members of the previous family is performed in 
Section 6. In Section 7 we describe the general solutions of a family of second-order Gambier 
equations in terms of particular solutions of other differential equations. We present a second 
quasi-Lie scheme for investigating second-order Gambier equations in Section 8, and conditions 
are given to be able to transform these equations into second-order Riccati equations. Those 
second-order Gambier equations that can be transformed into second-order Riccati equations 
are integrated in Section 9. Finally, Section 10 is devoted to summarising our main results. 
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2 Fundamentals 

Let us survey the fundamental results to be used throughout the work (see [15, 16, 17, 18] for 
details). In general, we hereafter assume all objects to be smooth, real, and globally defined on 
linear spaces. This simplifies our exposition and allows us to avoid tackling minor details. 

Given the projection ir : (t, x) G R x W 1 \— > x G R n and the tangent bundle projection 
r : TR n — > R n , a t-dependent vector field on R n is a mapping X : R x R n — > TR n such that 
r o X = ir. This condition entails that every t-dependent vector field X gives rise to a family 
{Xt}teR of vector fields Xt : x G 1" i— > X(t,x) G TR n and vice versa. We call minimal Lie 
algebra of X the smallest real Lie algebra V containing the vector fields {Xt}t^- Given a finite- 
dimensional M-linear space V of vector fields on R n , we write V(C°°(R)) for the C°° (M)-module 
of t-dependent vector fields taking values in V. 

An integral curve of X is a standard integral curve 7 : R — > R x R n of its suspension, i.e. 
the vector field X = d/dt + A(t, 2;) on E x M n . Note that the integral curves of X of the form 
7 : t G M — > (t, x(t)) elxl™ are the solutions of the system 



the referred to as associated system of X. Conversely, given such a system, we can define 
a t-dependent vector field on IR n [16] 



whose integral curves of the form (t,x(t)) are the solutions to (2.1). This justifies to write X 
for both a t-dependent vector field and its associated system. 

We call generalised flow a map g : (t,x) G R x R n h-» <7t(x) G M™ such that go = Id^n. 
Every t-dependent vector field X can be associated with a generalised flow g satisfying that the 
general solution of X can be written in the form x(t) = gt(xo) with xq G R n . Conversely, every 
generalised flow defines a vector field by means of the expression [15] 

ds s=t 

Generalised flows act on t-dependent vector fields [18]. More precisely, given a generalised 
flow g and a t-dependent vector field X, we can define a unique t-dependent vector field, 
g+X, whose associated system has general solution x(t) = g t (x(t)), where x(t) is the general 
solution of X. In other words, every g induces a t-dependent change of variables x(t) = g t (x(t)) 
transforming the system X into g+X. Indeed, g can be viewed as a diffeomorphism g : (t, x) G 
R x R n 1 — y (t,gt(x)) G R x R n , and it can easily be proved that g+X is the only t-dependent 
vector field such that g+X = g*X, where 5* is the standard action of the diffeomorphism g on 
vector fields (see [18]). 

Among all t-dependent vector fields, we henceforth focus on those whose associated systems 
are Lie systems. The characteristic property of Lie systems is to possess a superposition rule [7, 
18, 46]. A superposition rule for a system X on R n is a map <I> : (u^, . . . , u^ m y, k\, . . . , k n ) G 
(R n ) m x R n 1— > <3?(ii(i), . . . , Ur m y, fci, . . . , fc n ) G M n allowing us to write its general solution x(t) as 



dx % 
~dt 



X l {t,x) 



i = l,...,n, 



(2.1) 




8=1 



x(t) = $(#(i)(t),..., :E( m )(i); &i,...,fc n ) 



for a generic family of particular solutions X(i)(t), . . . , X( m )(t) and a set of constants k\, . . . , k. 
to be related to initial conditions. 



■n 
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The celebrated Lie-Scheffers theorem [46, Theorem 44] states that a system X possesses 
a superposition rule if and only if it is a i-dependent vector field taking values in a finite- 
dimensional real Lie algebra of vector fields, termed Vessiot-Guldberg Lie algebra [34, 53]. In 
other words, X is a Lie system if and only if V x is finite-dimensional [9]. This is indeed the 
main reason to define V x [14]. 

To illustrate the above notions, let us consider the Riccati equation [35] 

^ = b 1 (t) + b 2 (t)x + b 3 (t)x 2 , (2.2) 

where bi(t), b 2 (t), b 3 (t) are arbitrary functions of time. Its general solution, x(t), can be obtained 
from an expression [38, 56] 

x{t) = <£>(x {1) (t),x i2) (t),x i3) (t);k), 

where A: is a real number to be related to the initial conditions of every particular solution, 
£(!)(£), X( 2 )(t), £(3)(t) are three different particular solutions of (2.2) and $ : I 3 x 1 -> IR is 
given by 



"(i)(^(2) ~ "(3)) ~ ku (2) (u (3) - n (1) ) 
* «(l),«(J)i«(3) * = 7 S T7 S • 



That is, the Riccati equations admit a superposition rule. Therefore, from the Lie-Scheffers 
Theorem, we infer that the i-dependent vector field X associated to a Riccati equation is such 
that V x is finite-dimensional. Indeed, 

X = (b 1 (t) + b 2 (t)x + b 3 (t)x 2 )^. 

Taking into account that X\ = d/dx, X 2 = xd/dx, X 3 = x 2 d/dx span a finite-dimensional real 
Lie algebra V of vector fields and X t = bi(t)Xi + b 2 (t)X 2 + b 3 (t)X 3 , we obtain that {X t }teR C V 
and V x becomes a (finite-dimensional) Lie subalgebra of V. 

The Lie-Scheffers theorem shows that just some few first-order systems are Lie sys- 
tems [9, 36]. For instance, this theorem implies that all Lie systems on the real line are, up 
to a change of variables, a particular case of a linear or Riccati equation [55]. Therefore, many 
other important differential equations cannot be studied through Lie systems (see [9] for exam- 
ples of this). In order to treat non-Lie systems, new techniques generalising Lie systems need 
to be developed. We here focus on the theory of quasi-Lie schemes [15, 20]. 

Definition 2.1. Let W, V be finite-dimensional real vector spaces of vector fields on W 1 . We 
say that they form a quasi-Lie scheme S(W, V) if: 

• W is a vector subspace of V. 

• W is a Lie algebra of vector fields, i.e. [W, W] C W. 

• W normalises V, i.e. [W, V] C V. 

Associated to each quasi-Lie scheme, we have the C°°(IR)-modules W(C°°(M)) and V(C°°(M)) 
of t-dependent vector fields taking values in W and V, respectively. Now, from the Lie alge- 
bra W, we define the group Q(W) of generalised flows of t-dependent vector fields taking values 
in W, the so-called group of the scheme. The relevance of this group is due to the following 
theorem [15]. 

Theorem 2.1. Given a quasi-Lie scheme S(W,V), every generalised flow ofQ(W) acts trans- 
forming elements of V(C°°(W)) into members of V (C°° (M)) . 
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In other words, the elements of the group of a scheme provide t-dependent changes of variables 
that transform systems of V(C°°(M)) into systems of this family. Roughly speaking, we can un- 
derstand this group as a generalisation of the t-independent symmetry group of a system: apart 
from transforming the initial system into itself, the transformations of the group of a scheme 
also may transform the initial system into one of the "same type" . For instance, given a Lie 
system X associated with a Vessiot-Guldberg Lie algebra V, then S(V, V) becomes a quasi-Lie 
scheme. The group Q{V) allows us to transform X into a Lie system with a Vessiot-Guldberg 
Lie algebra V. This can be employed, for example, to transform Riccati equations into Riccati 
equations that can be easily integrated, giving rise to methods to integrate Riccati equations. 

In order to illustrate previous notions, we now turn to proving that quasi-Lie schemes allow 
us to cope with Abel equations of first-order and first kind [11], i.e. 

^ = h(t) + b 2 (t)x + b 3 (t)x 2 + b 4 (t)x 3 , (2.3) 

with bi(t), . . . , 64 (i) being arbitrary t-dependent functions. Indeed, if we fix W = (d/dx, xd/dx), 
it can be proved that SiW, V) and V = (d/dx, xd/dx, x 2 d/dx, x^d/dx) is a quasi-Lie scheme and 
X G V(C°°(R)) for every X related to an Abel equation (2.3). The elements of Q(W) transform 
Abel equations into Abel equations and geometrically recover the usual t-dependent changes 
of variables used to study these equations. This was employed in [11] to describe integrability 
properties of Abel equations. 

Given a quasi-Lie scheme S(W,V), certain systems in V(C°°(M)) can be mapped into Lie 
systems admitting a Vessiot-Guldberg Lie algebra contained in V. This enables us to study 
the transformed system through techniques from the theory of Lie systems and, undoing the 
performed transformation, to obtain properties of the initial system under study [15]. 

Definition 2.2. Let S(W, V) be a quasi-Lie scheme and X a t-dependent vector field in 
V(C°°(M)), we say that A is a quasi-Lie system with respect to S(W,V) if there exists a gene- 
ralised flow g E Q(W) and a Lie algebra of vector fields Vq C V such that g+X G Vo(C°°(M)). 

3 A new quasi-Lie scheme 

for investigating second-order Gambier equations 

The Gambier equation [30, 32] can be described as the coupling of two Riccati equations in 
cascade, which can be given in the following form 

dy 2 

— = - y j r a 1 y + a 2 , 
dx 2 

— = aox + nyx + a, 
dt 

where n is an integer, a is a constant, which can be scaled to 1 unless it happens to be 0, 
and ao, oi, a 2 are certain functions depending on time. The precise form of the coefficients 
of the Gambier equation is determined by singularity analysis, which leads to some constraints 
on ao, oi and a 2 [30]. For simplicity, we hereafter assume ao(0) ^ 0. Nevertheless, all our results 
can easily be generalised for the case ao(0) = 0. 

If n 7^ 0, we can eliminate y between the two equations above, which gives rise to the referred 
to as second-order Gambier equation [32, 39, 40, 44], i.e. 

d 2 x n — 1 f dx \ 2 (n + 2) dx dx (n — 2) dx 
1 - +a - -x— + a 1 — -ct- v 



dt 2 xn \dt J n dt dt nx dt 

aZ o / dao \ o ( a \ 0-2 

x + — aoGq x + a 2 n — 2a$— x — a±a . (3.1) 

n \ dt ) \ nJ nx 
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The importance of second-order Gambier equations is due to their relations to remarkable 
differential equations such as second-order Riccati equations [31, 35], second-order Kummer- 
Schwarz equations [5, 16] and Milne-Pinney equations [32]. Additionally, by making appropriate 
limits in their coefficients, Gambier equations describe all the linearisable equations of the 
Painleve-Gambier list [32] . Several particular cases of these equations have also been studied in 
order to analyse discrete systems [40]. 

Particular instances of (3.1) have already been investigated through the theory of Lie systems 
and quasi-Lie schemes. For instance, by fixing n = —2, a = a\ = 0, and ao to be a constant, the 
second-order Gambier equation (3.1) becomes a second-order Kummer-Schwarz equation (KS2 
equation) [5, 16] 

d 2 x 3 f dx\ 2 o . , . 

W = 2x[dt) -W + M*K (3-2) 

where we have written cq = — Og/4, with Co a non-positive constant, and u(t) = — 02(f) so as 
to keep, for simplicity in following procedures, the same notion as used in the literature, e.g. 
in [16]. The interest of KS2 equations is due to their relations to other differential equations 
of physical and mathematical interest [16, 24, 32]. For instance, for x > the change of 
variables y = 1/y/x transforms KS2 equations into Milne-Pinney equations, which frequently 
occur in cosmology [32]. Meanwhile, the non-local transformation dy/dt = x maps KS2 equations 
into a particular type of third-order Kummer-Schwarz equations, which are closely related to 
Schwarzian derivatives [4, 16, 43]. Additionally, KS2 equations can be related, through the 
addition of the new variable dx/dt = v, to a Lie system associated to a Vessiot-Guldberg Lie 
algebra isomorphic to st(2,R), which gave rise to various methods to study its properties and 
related problems [16]. 

If we now assume n = 1 and a = in (3.1), it results 

d 2 x . dx 23 f da-o \ 

— = (ai + 3a x) — - a x + I — - a ai I x + a 2 x, 

which is a particular case of second-order Riccati equations [23, 31] that has been treated through 
the theory of quasi-Lie schemes and Lie systems in several works [10, 13, 29]. Furthermore, 
equations of this type have been broadly investigated because of its appearance in the study of 
the Backhand transformations for PDEs, their relation to physical problems, and the interest of 
the algebraic structure of their Lie symmetries [8, 23, 31, 37, 42]. 

In view of the previous results, it is natural to wonder which kind of second-order Gambier 
equations can be studied through the theory of quasi-Lie schemes. To this end, let us build up 
a quasi-Lie scheme for analysing these equations. 

As usual, the introduction of the new variable v = dx/dt enables us to relate the second-order 
Gambier equation (3.1) to the first-order system 

dx 

dv in - 1) v 2 (n + 2) (n - 2) v al % 

— = h a - -xv + aiv - a- ^-x 6 

dt n x n n x n 

( dao \ q / _ a\ a 2 

+ — a ai ) x + a 2 n - 2a - x - a±a , 

\ at J V nJ nx 

which is associated to the t-dependent vector field on TMo ; with Mq = K \ {0}, given by 



ox 



(n - 1) v 2 (n + 2) (n-2)v al n 

h a - -xv + aiv - a- -x 6 

n x n n x n 



A Quasi-Lie Schemes Approach to Second-Order Gambier Equations 



7 



dag \ 2 , / o a \ 0-2 



+ — aoa\ ) x + ( a2« — 2ao— x — ai<7 

at /V n/ ra 



St;' 



termed henceforth the Gambier vector field. To obtain a quasi-Lie scheme for studying the above 
equations, we need to find a finite-dimensional M- linear space Vg such that X G Vg^C 00 ^)) 
for all do, ai, a2, f and n. Observe that X can be cast in the form 

X = ^ b a (a , -"p ai, a 2 , cr, n j F a , (3.3) 

with 6i, . . . , 610 being certain t-dependent functions whose form depends on the functions ao, 
dao/dt, ai, a 2 and the constants cr and n. More specifically, these functions read 

h = l, b 2 = , h = a , 64 = 01, h = -a , 

n n n 

c>6 = , 07 = -— - a ai, 6 8 = a2W - 2a -, 69 = -aia, 610 = 

n dt n n 

and 

v 8 v y2 d v 9 v 9 v v 9 

Y 1 = v—, Y 2 = — — , Y 3 = xv — , F 4 = u— , F 5 = — — , 

ax x ov ov x ov 

v 3 9 2 ^ ^ d a 1 3 

OV OV OV OV X ov 

For convenience, we further define the vector field 
ox 

which, although does not appear in the decomposition (3.3), will shortly become useful so as to 
describe the properties of Gambier vector fields. 

In view of (3.3), it easily follows that we can choose Vg to be the space spanned by Yi, . . . , Yu. 
It is interesting to note that the linear space Vq is not a Lie algebra as [I3 , Yq] does not belong 
to Vq- Moreover, as 

j'-times 

^y 6 = \y 3 ,[y 3 X..,[Y3,y 6 ] . . .}}} = (-i)V+ 3 ^, j e N, 

3 OV 

there is no finite-dimensional real Lie algebra V D Vg such that X G V(C°°(M)). Hence, X is 
not in general a Lie system, which suggests us to use quasi-Lie schemes to investigate it. 

To determine a quasi-Lie scheme involving Vg, we must find a real finite-dimensional Lie 
algebra Wq C Vg such that [Wg, Vg] C Vg- In view of Table 1, we can do so by setting 
Wq = {Y4,Yg,Yu) , which is a solvable three-dimensional Lie algebra. In fact, 

[Y4, Y 8 ] = -Y 8 , [Y 4 , Y n ] = 0, [Y 8 , Y n ] = -Y 8 . 

In other words, we have proved the following proposition providing a new quasi-Lie scheme to 
study Gambier vector fields and, as shown posteriorly, second-order Gambier equations. 

Proposition 3.1. The spaces Vg = (Y\, . . . ,Yn) and Wg = (Ya,Y%,Yh) form a quasi-Lie 
scheme S(Wg, Vg) such that X £ Vg(C°°(J&)) for every Gambier vector field X. 

Recall that Yn is not necessary so that every Gambier vector field takes values in Vg- Hence, 
why it is convenient to add it to Vg? One reason can be found in Table 1. If Vg had not 
contained Yn, then S(Wg, Vg) would have not been a quasi-Lie scheme as Wg ^ Vg- In 
addition, Yg could not belong to Wg neither, as [Yg, Yl] G Vg provided Yn — Y4 G Vq. Hence, 
including Yn in Vq allows us to choose a larger Wg C Vq- In turn this gives rise to a larger 
group Q(Wg), which will be of great use in following sections. 
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Table 1. Lie brackets \Y U Yj] with i = 4,8, 11 and j = 1, ... ,11. 



[v] 


Yt 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 


Y 7 


Y s 


Yg 


Y 10 


Y n 


Y 4 


Yt 


Y 2 











-Y 6 


-Y 7 


-Ys 


-Y 9 


-Y w 





Y s 


Y n ~ Y 4 


2Y 4 


Y 7 


Y 8 


^9 

















-Ys 


Y n 


-Y x 


-Y 2 


Y 3 





-Y 5 


3^6 


2Y 7 


Y 8 





-Y 10 






4 Transformation properties of second-order Gambier equations 

Remind that Theorem 2.1 states that the i-dependent changes of variables associated to the 
elements of the group G(W) of a quasi-Lie scheme S(W, V) establish bijections among the t- 
dependent vector fields taking values in V. This may be of great use so as to transform their 
associated systems into simplified forms, e.g. in the case of Abel equations [11]. We next show 
how this can be done for studying the transformation of second-order Gambier equations into 
simpler ones whose corresponding coefficients a\ vanish [30]. This retrives known results from 
a geometrical viewpoint and shows that certain second-order Gambier equations cannot be 
transformed into simpler ones, solving a small gap performed in [30]. 

As the vector fields in Wq span a finite-dimensional real Lie algebra of vector fields on TRoj 
there exists a local Lie group action (p : G x TMo ~~ ^ TIRo whose fundamental vector fields are 
the elements of Wg- By integrating the vector fields of Wq (see [12] for details), the action can 
easily be written as 



y " , -H«)) = U + 8v)> Wh6re 9eTd '{(l °S 



a, 5 G E + , 7 G 



The theory of Lie systems [9, 17] states that the solutions of a system associated to a t-dependent 
vector field taking values in the real Lie algebra Wq are of the form (x(t),v(t)) = (p(h(t), (xq,vq)), 
with h{t) being a curve in with /i(0) = e. Therefore, every g G Q(Wg) can be written as 
gt(-) = (p(h(t), •) for a certain curve h(t) in G with h(0) = e. Conversely, given a curve h(t) in 
Td with h(0) = e, the curve (x(t),v(t)) = cp(h(t), (xq,vq)) is the general solution of a system of 
W G (C°°(R)), which leads to a generalised flow g t {-) = <p{h{t), •) of G{W G ) [9, 17]. Hence, the 
elements of Q{Wq) are generalised flows of the form 

g h V(t,x,v) = <p(h(t),x,v), 

for h(t) any curve in Td with h(0) = e. Observe that every h(t) is a matrix of the form 



hit) 



a(t) 
7 (t) 5(t) 



where a(t), 5(t) and j(t) are t-dependent functions such that a(0) = 5(0) = 1 and 7(0) = 0, 
because h(0) = e, and a(t) > 0, S(t) > as h(t) G for every i 6 1. Hence, every element 
of Q(Wg) is of the form 

g <*(t)a(t)m(t,x,v) = g h ®(t,x,v) = (t,a(t)x,j(t)x + S{t)v). (4.1) 

Theorem 2.1 implies that for every g G Q{Wg) an d Gambier vector field X G Vg(C°° (M.)), 
we have g+X G Vg(C°°(M)). More specifically, a long but straightforward calculation shows 
that 

11 

g+X = Y J b*(t)Y a , (4.2) 

a=l 
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where the functions b a = b a (t), with a = 1, . . . , 11, are 

- a - n — 1 a - anin + 2) - 1 d5 (2 — 71)7 
h = j, b 2 = -, 6 3 = — b i = a 1 + -—+ y - 

no na at no 

T (2-n)aa T aZS T 1 / ^da n + 2 

h = , h = Sr, h = —~ 5— - a aid a 7 

n na° ar \ at n 

- 5 ( 2a 7 7 2 7 dS 1 drf 
8 a \ 2 n° 5 1 n5 2 5 2 dt 5 dt 

- ( . (2 — 71)7 \ - a 2 a5 - Ida 7 fJn . 
b 9 = -a[a 1 5+± ^ , 610 = , b n = -- j 7-j- (4.3) 

\ n J n a dt 

Let us use this so as to transform the initial Gambier vector field into another one that is related, 
up to a t-reparametrisation r = r(t), to a second-order Gambier equation with r-dependent 
coefficients olq, di, d~2, a constant a and an integer number n. Additionally, we impose d% = 
in order to reduce our initial first-order Gambier equation into a simpler one. In this way, we 
have 



g*X = m Yl + ^Y 2 + d ^Y 3 -a—Y 5 
n n n 

On dan / a\ a 2 

4y 6 + — °y 7 + (a 2 n - 2d Q A y 8 - —y 10 



(4.4) 



for a certain non- vanishing function £(t) = dr/dt. Therefore, 64 = 69 = b\\ = 0, i.e. 



Id5 (2-n) 7 . , 

aLs+^A=0, (4.6) 

(4.7) 

a at 

As we want our method to work for all values of a, e.g. <j/0, equation (4.6) implies 

M+ (2-n) 7 =Q _ 
n 

As 5 > 0, the above equation involves that a\ = for n = 2. In other words, we cannot 
transform a Gambier vector field with a\ 7^ into a new one with d\ = through our methods 
if ri = 2 and a 7^ 0. In view of this, let us assume that n 7^ 2. 

From (4.5) and (4.8), and using again that 5 > 0, we infer that dS/dt = 0. As (5(0) = 1, then 
5 = 1. Plugging the value of 5 into (4.7) and (4.8), we obtain 

1 da na\ ( f l na\ A na± 



a = exp (i i^h dt ) ' 7 



a (it n — 2 \J n — 2 J n — 2 

which fixes the form of g mapping a system X into (4.4). Bearing previous results in mind, we 
see that the non-vanishing i-dependent coefficients (4.3) become 

r 7 n ~ 1 7 oo (n + 2) - (2 - n)o-a - 

61 = a, 6 2 = a, 03 = , 65 = , b 6 = a , 

n rta n na- 3 

- d /ao\ - 1/ 2<rao n(n — 1) 2 n dai\ - a 2 a 
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Comparing this and (4.4), we see that to transform the initial first-order Gambier equation into 
a new one through a i-dependent change of variables (4.1) and a i-reparametrisation dr = adt 
requires £ = a. The resulting system reads 

dx 

dv = (n-1) * + ^ (n + 2)^ _ ^(n-2) v _ ^ 

dr n x a l n n x a^n 

2^ 



, d(a /a j , , 

H : -x + a^n — iaoa — x 



dr V nJnx 



if (n — 1) 2 1 dai 



where 

Therefore, redefining n = n, ao = ag/a 2 and cj = u, we obtain that the above system is 
associated to a second-order Gambier equation with d\ = 0. Meanwhile, as g induces a t- 
dependent change of variables given by 

da 

x = ax, v = —x + v, 
dr 

we see that this i-dependent change of variables can be viewed as a t-dependent change of 
variables x = ax along with a t-reparametrisation t — > t such that v = dx/dr. Indeed, 

dx da dx da 
x = ax — = —x + a— = —x + v = v. 

dr dr dr dr 

Furthermore, these transformations map the initial second-order Gambier equation with 

into a new one with d\ = 0, i.e. depending only on two functions ao and 02 and the constants a 

and n. We can therefore formulate the following result. 

Proposition 4.1. Every second-order Gambier equation (3.1) with n ^ 2 can be transformed via 
a t-dependent change of variables x = a{t)x and a t-reparametrisation t = r(t), with dr = adt, 
a(0) = 1 and 

1 da n 
a dt n — 2 

into a second-order Gambier equation whose a\ vanishes and n, a remain the same. 

In the above proposition, we excluded second-order Gambier equations with n = 2 as we 
noticed that in this case the proof of this proposition does not hold: we cannot transform an 
initial Gambier vector field with a 7^ and a\ ^ into a new one with d\ = 0. Moreover, it 
is easy to see that the transformation provided in Proposition 4.1 does not exist for n = 2 and 
ai 7^ 0. This was not noticed in [30], where this transformation is wrongly claimed to transform 
any Gambier equation into a simpler one with a\ = 0. In view of this, we cannot neither ensure, 
as claimed in [30], that second-order Gambier equations are not given in their simplest form. 



5 Quasi-Lie systems and Gambier equations 

The theory of quasi-Lie schemes mainly provides information about quasi-Lie systems, which 
can be mapped to Lie systems through one of the transformations of the group of a quasi-Lie 
scheme. This allows us to employ the techniques of the theory of Lie systems so as to study the 
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obtained Lie systems, and, undoing the performed i-dependent change of variables, to describe 
properties of the initial system [15]. 

Motivated by the above, we determine and study the Gambier vector fields X G V(C°°(]R)) 
which are quasi-Lie systems relative to S(Wq, Vg)- This task relies in finding triples (g, X, Vq), 
with g G Q(W G ) and Vq being a real Lie algebra included in Vq in such a way that g+ X G 
V (C°°(R)). 

One of the key points to determine quasi-Lie systems is to find a Lie algebra Vo- In the 
case of Gambier vector fields, this can readily be obtained by recalling that certain instances of 
second-order Gambier equations, e.g. (3.2), are particular cases of KS2 equations. By adding 
a new variable v = dx/dt to (3.2), the resulting first-order system becomes a Lie system (see [16]) 
related to a three-dimensional Vessiot-Guldberg Lie algebra Vo C Vq of vector fields on TMo 
spanned by 

X 1 = 2x—, X 2 = x— + 2v—, X 3 = v— + (---2c x 3 )—, (5.1) 



dv' ' dx dv : ' dx \2 x J dv 

i.e. 

X 1 = 2Y 8 , X 2 = Y ll + 2Y A , X 3 =Y 1 + ~Y 2 -2c Y fi . 

Consequently, it makes sense to look for Gambier vector fields X G V(C°°(M)) that can be 
transformed, via an element g G G(Wg), into a Lie system g+X G Vb(C°°(R)), i.e. 

g+ X = 2f(t)Y 8 + g(t)(Y n + 2Y 4 ) + h(t) (y x + ^Y 2 - 2c Y 6 

for certain t-dependent functions /, g and h. Comparing the expression of g+X given by (4.3) 
and the above, we find that g+X G Vb(C°°(lR)) if and only if 

h = h = h = b 9 = b l0 = (5.2) 

and 

_ 3_ _ 

6 4 = 26n, b 2 = -61, b 6 = -2c 6i. (5.3) 

From expressions (4.3) and remembering that a > and 5 > 0, we see that condition 610 = 
implies that a = 0. This involves, along with (4.3), that 65 = 69 = 0. Meanwhile, from 
b 2 = 3&i/2 we obtain n = —2, which in turn ensures that 63 = 0. Bearing all this in mind, 
67 = reads 

dan 

a ai - -J- = 0. (5.4) 

Above results impose restrictions on the form of the Gambier vector field X to be able to be 
transformed into a Lie system possessing a Vessiot-Guldberg Lie algebra Vo via an element 
g G G(Wg)- Let us show that the remaining conditions in (5.2) and (5.3) merely characterise 
the form of the i-dependent change of variables g and the coefficient cq appearing in (3.2). 
Conditions 2b\i = &4 and be = — 2co&i read 

r 2 4 

d , a 2 a . . 

^!ogy = ai, a = -4c -^-. (5.5) 

Using the first condition above, the relation (5.4), and taking into account that ck(0) = 5(0) = 1 
and 5, a > 0, we see that 



^(a 2 <5 _1 ) _ 1 da a 2 _ a 



a 2 5 1 ao dt 5 gso(0) 
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Using the second equality in (5.5), we obtain 
4c = -a (0) 2 , 

which fixes cq in terms of a coefficient of the initial second-order Gambier equation. 

Concerning the i-dependent coefficients b\, . . . , b\ii the non- vanishing ones under above con- 
ditions, i.e. b\, h, t>4, &6; bg and &n, can readily be obtained through relations (5.3) and 

t _« t _^_f_r > l_dao _ j^d5_ ld_7_\ 7 Ida 7 

1 ~ 5 ' 8 ~ a V Ja (it 2<5 2 <5 2 dt + 5 dt J ' 11 " a dt 5" 

In other words, we have proved the following proposition. 

Proposition 5.1. A Gambier vector field X is a quasi-Lie system relative to S(Wg, Vg) that 
can be transformed into a Lie system taking values in a Lie algebra of the form Vq if and only 
if a$ai = dao/dt, n = —2 and a = 0. Under these conditions, the constant cq appearing in X3 
becomes Co = — ao(0) 2 /4. Additionally, a transformation g(a(t), r y(t),8(t)) G G(Wg) maps X 
into g+X G Vb(C°°(M)) if and only if a 2 /5 = ao/ao(0). More specifically, g+X reads 

a f Ida 7\ 5 ( 7 da 7 2 7 <W , 1 d-y\ , . 

9^ = -X 3+ ^----jX 2 + -[-2a 2 -—— + W2 - T2 - + --jX 1 . (5.6) 

The above proposition allows us to determine the transformations g £ Q(Wq) ensuring that 
a Gambier vector field and its related second-order Gambier equation can be mapped, maybe 
up to a t-reparametrisation, into a new Gambier vector field related to a KS2 equation. Indeed, 
from (3.2), we see that to do so, we need to impose 

g+X = H(t)(X 3 + cj(t)X 1 ), 

for a certain function ui(t) and a function such that 7^ for all t EK. Comparing (5.6) 
with the above expression, we see that 

a at 

which, in view of the fact that <7 a (*)>TW> 5 (*) satisfies 
ao{0) 

enables us to determine the searched transformations. In fact, fixed a non- vanishing ^-dependent 
function a with a(0) = 1, the above conditions determine the values of 5 and 7 of g. 
Now, a t-dependent reparametrisation 

* a 

T dt' (5.9) 


transforms the system associated to g+X into a new system related to X3 + (jj(t(r))X\, where 
w is given by 

W = ~2c7 [ 2a2 + aoS^t " W ~ STt + P dt ) • (5 - 10) 

Proposition 5.2. Every Gambier vector field X satisfying ao°i = dao/dt, n = —2 and <r = 
is a quasi-Lie system relative to S(Wg, Vg) that can be transformed into a Lie system associa- 
ted to X3 + u)(t)Xi, with Co = — ao(0) 2 /4 and u(t) given by (5.10), through a transformation 
^a(*)i7(*)/(t) £ G(Wg), whose coefficients are given by any solution of (5.7) and (5.8), and the 
t-dependent reparametrisation (5.9). 
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Note that the transformation g Q (*)'TM> 5 M ca n be viewed as a t-dependent change of variables 
x = ax and a t-reparametrisation dr = a/Sdt, i.e. 

x = ax, v = 7X + 5v, 

and, in view of the first condition in (5.7), we see that 

dx da dx 

— = —x + a— = 7X + ov = v. 

dr dr dr 

From this and Proposition 5.2, we obtain the following proposition about second-order Gambier 
equations. 

Corollary 5.1. Every second-order Gambier equation of the form 

d 2 x 3 ( dx\ 2 1 dan dx a 2 , n ,„ 

' H j^-r + — x ~ 2a 2X (5.11) 



dt 2 2x V dt J a dt dt 2 
can be mapped into a KS2 equation 

d 2 x 3 (dx\ 2 1 2 _ 3 5 2 ( 7 dan lj 2 1 dj 7 d5\ . 

+ -a (0) x° ? 2a 2 H j— - - 7-77 + 7^-7- x (5.12) 



dr 2 2x \dr J 2 w a 2 V a <5 2 J 2 5 (5 2 dt 

through a t-dependent change of variables x(t) = a{t)x{t), where a{t) is any positive function 
with a(0) = 1, and a t-reparametrisation r(t) with dr = ajbdt, with 5 and 7 being determined 
from a by the relations 

2 a (0) aa (0)da 

d = a , 7= — . (5.13) 

ao an dt 

6 Constants of motion for second-order Gambier equations 

In this section, we obtain constants of motion for second-order Gambier equations. We do so 
by analysing the existence of t-independent constants of motion for systems g+ X £ Vb(C°°(]R)), 
with X being a Gambier vector field. By undoing the t-dependent change of variables g, this 
leads to determining constants of motion for a Gambier vector field X and its corresponding 
second-order Gambier equation. 

Let F : TIRo — )■ M be a t-independent constant of motion for g+X, we have (g+X)tF = 
for all f G K. This involves that F is a t-independent constant of motion for all successive Lie 
brackets of elements of {(g+X)t}teM. as wen as their linear combinations. In other words, F is 
a common first- integral for the vector fields V 9 * x C Vn. 

When u(t) is not a constant, it can be verified that V 9 * x = Vn. Thus, if F is a first-integral 
for all these vector fields, it is so for all vector fields contained in the generalised distribution 
V, p = ((Xi) p , (X2) p , (A3) p ), with p £ TIRo- Hence, dF p G V°, i.e. dF p is incident to all vectors 
of T>p. In this case, T> p = T p TlRo for a generic p G TRo> which implies that dF p = at 
almost every point. Since F is assumed to be differentiable, we have that F is constant on each 
connected component of TMo and g+X has no non-trivial t-independent constant of motion. 

If u(t) = A for a real constant A, then dim2? p = 1 at a generic point and it makes sense to 
look for non-trivial t-dependent constants of motion for g+X. In view of (5.10), this condition 
implies 



a (0) 2 a 2 



2a 2 



1 dan d log a 1 / d log a \ 2 d 2 log a 
° 2 + ao^t dt 2 V dt ) dt 2 
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The function F can therefore be obtained by integrating 
dF [3v 2 a (0) 2 » , \5F 

In employing the characteristics method [10], we find that F must be constant along the integral 
curves of the so-called characteristic system, namely 

dx dv 



Let us focus on the region with v ^ 0, i.e. O = {(x,v) G TIRo I v 7^ 0}. We obtain from the 
previous equations that 

dv 3 v an(0) 2 x 3 , a; 

— = — + + 2A-. 

dx 2 x 2 v v 

Let us focus on the case x > 0; the other case can be obtained in a similar way and leads to the 
same result. Multiplying on right and left by v/x and defining w = v 2 and z = x 2 , we obtain 

dw _3w a (0) 2 

-3— — 1 ~ Z + ZA. 

dz 2 z 2 

As this equation is linear, its general solution can be easily derived to obtain 

W ( Z ) = (ao(O)V/ 2 - 4AZ- 1 / 2 + £)z 3 / 2 , 

for an arbitrary real constant £. From here, it results 

^2 v 2 4A t 
-a (0 yx + ^A = £. 

x a x 

Consequently, F is any function of the form F = F(£), for instance, 

v 2 4A 

F = -a (0) 2 x + ^ + — , (x,v)€0. 

X a X 

In principle, F was defined only on O. Nevertheless, as this region is an open and dense subset 
of TIRo) an d in view of the expression for F, we can extend it differentiably to TIRo m a unique 
way. Since F is a constant of motion on O, it trivially becomes so on the whole TIRo- 
Summarising, we have proved the following. 

Theorem 6.1. A second-order Gambier equation 

d 2 x 3 f dx\ 2 1 dan dx a 2 , n . 

■7^ = - — + — + -^x 3 - 2a 2 x, 6.1 

dt 2 2x \ dt J a dt dt 2 ' y ' 

admits a constant of motion of the form 

F = -a (0) 2 x + ^ + —, (6.2) 
x- 3 x 

where X is a real constant, x = ax and v = 5dx/dt + jx, with a being a particular positive 
solution with a(0) = 1 of 

d 2 loga 2Aan „ 1 dan d log a 1 /<iloga\ 2 ,„ „. 

+ 2a 2 + -if— yf -° (6.3) 



di 2 a (0) 2 a 2 a dt dt 2\ dt 

and 7, 5 are determined from a by means of the conditions (5.13). 
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The above proposition can be employed to derive a constant of motion for certain families 
of second-order Gambier equations. For instance, if we start by a Gambier equation (6.1) with 
a2 = — Aag/ao(0) 2 , for a certain real constant A, then a = 1 is a solution of (6.3). In view 
of (5.13), 7 = and 5 = ao(0)/ao. Therefore, Theorem 6.1 establishes that the second-order 
Gambier equation (6.1) admits a constant of motion 

, ,i oo(0) 2 / 'dx\ 2 4A . 
F = -o (0) 2 x + ^-^ — +— . (6.4) 



clqX 3 \dt J x 

Consider now a general second-order Gambier equation (6.1) and let us search for a constant 
of motion (6.2) with A = 0. In this case, (6.3) can be brought into a Riccati equation 

dw 1 dao 1 2 

— = 2a 2 + —w - -w , 

at ao dt 2 

where w = d log a/dt, whose solutions can be investigated through many methods [19, 22]. The 
derivation of a particular solution provides a constant of motion for the second-order Gambier 
equation (6.1) that can be obtained through the previous theorem. Additionally, this particular 
solution can be used to obtain the general solution of the Riccati equation [9, 19], which in turns 
could be used to derive new constants of motion for the second-order Gambier equation (6.1). 

Note that all the above procedure depends deeply on the fact that A is a constant. Recall 
that if A is not a constant, then g+X does not admit any ^independent constant of motion. 
Nevertheless, other methods can potentially be applied in this case. For instance, using that 
S(Vo,Vo) is a quasi-Lie scheme, we can derive the group Q(Vq) of this scheme and to use an 
element h £ Q{Vq) to transform g+X into other Lie system h+g+X of the same type, e.g. 
another of the form h+g+X = £,2{t){c\Xi + C2A2 + C3X3), with ci, 02,03 G M, the vector fields 
X\, X2, A3 are those of (5.1), and £2^) is any t-dependent nonvanishing function. As this 
system is, up to a t-parametrization, a t-independent vector field, it admits a local t-independent 
constant of motion. By inverting the t-dependent changes of variables h and g, it gives rise to 
a t-dependent constant of motion of X and the corresponding second-order Gambier equation. 



7 Second-order Gambier and Milne— Pinney equations 

Consider the KS2 equation (5.12) with x > (we can proceed analogously for the case x < 0). 
The change of variables x = 1/y 2 , with y > 0, transforms it into a Milne-Pinney equation [52] 

d 2 y . . «o(0) 2 r71 s 

— 2 =-u{t{r))y-^. (7.1) 

These equations admit a description in terms of a Lie system related to a Vessiot-Guldberg Lie 
algebra isomorphic to sl(2,M) [12, 56]. This was employed in several works to describe their 
general solutions in terms of particular solutions of the same or others differential equations, 
e.g. Riccati equations and t-dependent harmonic oscillators [9]. 

Previous results allow us to describe the general solution of (5.11) in terms of particular 
solutions of Riccati equations or t-dependent frequency harmonic oscillators. Indeed, in view of 
Corollary 5.1, these equations can be transformed into a KS2 equation through a t-dependent 
change of variables x(t) = a(t)x(t) and a t-reparametrisation dr = (a/5) dt. In turn, y = l/\/S 
transforms (5.12) into (7.1), whose general solution y(t) can be written as [12] 



y(r) = ^hz 2 (T) + k 2 z 2 2 {T) + 2C(h, fe 2 , W) Zi (t)z 2 (t), 
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where C 2 (fci, k 2 , W) = k%k 2 + ao(0) 2 / (4W 2 ), the functions Zi(r), z 2 {r) are two linearly indepen- 
dent solutions of the system 

d 2 z 

and W is the Wronskian related to such solutions. Inverting previous changes of variables, the 
general solution for any second-order Gambier equation (6.1) reads 

x(t) = cT 1 [fciz?(r(t)) + k 2 z 2 2 (T(t)) ± 2^1*2 + ao(0)y(4W 2 )z 1 (T(t))z 2 (T(t))' ] 
Therefore, we have proved the following proposition: 

Proposition 7.1. The general solution of a second-order Gambier equation 
d 2 x 3 ( dx\ 2 1 dao dx a 2 , 3 



, , + r^-p + —x* - 2a 2 x, (7.2) 

dt 2 2x \dt J a dt dt 2 ' y ' 

can be brought into the form 



x(t) = a' 1 [fci*?(r(*)) + k 2 z 2 2 (T(t)) ± 2y/hk 2 + a (0) 2 / 1 {AW 2 ) Zl {T(t))z 2 {T(t))\ , 
where z\(t) and z 2 {t) are particular solutions of the r-dependent harmonic oscillator 

d 2 z 5 2 ( 7 da l7 2 ldj -fdS\ 

— 2 = -u(t(r))z = -— 2 [2a 2 + — — -__--_ + __j z> 

and W = z\ dz 2 /dr — z 2 dz\jdr, with dr = a/5 dt and a, 5 and 7 certain t-dependent functions 
satisfying (5.13). 

Many other similar results can be obtained in an analogous manner. For instance, the theo- 
ry of Lie systems was also used in [12] to prove that the general solution of a Milne-Pinney 
equation (7.1) can be written as 



y(j) 



[fci(xi(r) - x 2 {r)) - k 2 ( Xl (r) - x 3 (r))] 2 - a (0) 2 [x 2 (r) - x 3 (r)]/4 



I (k 2 - ki)(x 2 {r) - x 3 (t))(x 2 (t) - x 1 {t))(x 1 (t) - x 3 (t)) 

where x±(t), x 2 (t) and x 3 (t) are three different particular solutions of the Riccati equation 

~ = - x 2 . 

Proceeding as before, we can describe the general solution of a second-order Gambier equa- 
tion (7.2) in terms of solutions of these Riccati equations. In addition, by applying the theory 
of Lie systems [19] to Milne-Pinney equations (7.1), we can obtain many other results about 
subfamilies of second-order Gambier equations. In addition, the relation between second-order 
Gambier equations and Milne-Pinney equations enables us to obtain several other results in 
a simple way. 

Proposition 7.2. The second-order Gambier equation (5.11) with 2a 2 = a 2 , > 0, i.e. 

d 2 x 3 / dx\ 2 1 dandx a 2 , q , ,_ . 

— ^ = — — + -H— + -^-x 3 -a 2 ) x, 7.3 

dt 2 2x \ dt J a dt dt 2 ' v ; 

can be transformed into the integrable Milne-Pinney equation 

d 2 y y 1 , 

— — — H = (7 4) 

dr 2 2 + 4y3 U ^ 

under the transformation x = 1/y 2 and a t-reparametrisation dr = a$dt. 
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Note 7.1. For second-order Gambier equations (5.11) with oq = 0, the transformations given 
in the above proposition map these equations into harmonic oscillators. 

Proposition 7.3. A constant of motion of (7.4) is given by 

2 



Imp 



dy 
dr 



2 4y 2 



Using y = x l l 2 transformation, we obtain a constant of motion for the special second-order 
Gambier equation (7.3) of the form 



he 



1 1 



4 x 3 al \dt J 



dx\' 



1 x 
2x~ + 4 



(7.5) 



Thus, we say the special equation (7.3), which yields a constant of motion, is an integrable 
deformation of the Milne-Pinney equation. 

Note 7.2. Recall that Theorem 6.1 can be applied to those equations (5.11) where 02 = 
— Aog/a^O) giving rise to the constant of motion (6.4). Observe that (7.3) is of this type 
with A = — ao(0) 2 /2. Then, the constant of motion (6.4) is, up to a multiplicative constant, the 
same as (7.5). Despite this, the above illustrates how certain results can readily be obtained 
through Milne-Pinney equations. 



8 A second quasi-Lie schemes approach 
to second-order Gambier equations 

Apart from our first approach, we can provide a second quasi-Lie scheme to study second-order 
Gambier equations. This one is motivated by the fact that some cases of these equations, namely 

d 2 x ( dx o\ „. . . . , , . / 9 dx\ 

W = - ( 3x dt + x ) + m + 9{t)x + h{t) [ x + -dt) ' 

for arbitrary f-dependent functions /, g and h, form a particular subclass of second-order Riccati 
equations that are SODE Lie systems [13, 16]. That is, when written as first-order systems 

dx 

Tt =V ' 

^ = - {Sxv + x 3 ) + f(t) + g{t)x + h(t) (x 2 + v) , (8.1) 

they become Lie systems. Indeed, these systems possess a Vessiot-Guldberg Lie algebra Vq 
isomorphic to sl(3, R) given by [10, 29] 



X 6 = 2x(v + x 2 )^- + 2(v 2 -x 4 )!-, x 7 = A_aJL x 8 = 2x^- + Av^-. 

ox ov ox ov ox ov 

In terms of these vector fields, (8.1) is the associated system to the i-dependent vector field 
X x + f(t)X 2 + ^g(t)(X 3 + X 7 ) + ^(-2X4 + X s ). 
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Table 2. Lie brackets [Yi, Yj] with i = 4, 8, 11 and j = 12, . . . , 17. 



[v] 


^12 


Yl 3 


Y u 


Y 15 


Y 16 


n 7 










Y u 





-Y 16 


*17 


^8 


-Yg 


-Y 7 


Y 13 - Y 3 


-Y 6 





213 




-Y12 


Y 13 





2Y 15 


4^16 






Therefore, it is natural to analyse which Gambier vector fields X can be transformed through 
a quasi-Lie scheme into one of these Lie systems. 

Despite the interest of the previous idea, the quasi-Lie scheme provided in the previous section 
cannot be employed to determine all quasi-Lie systems of this type. The reason is that Vo is 
not included in Vq and, therefore, there exists no g S Q(Wg) such that g+X G Vb(C°°(M)) C 
Vg(C°°(M)). To solve this drawback, we now determine a new quasi-Lie scheme S(Wg, Vq) such 
that the space V' G contains Vq + Vq. This can be done by defining V' G as the M- linear space 
generated by the elements of Vq and the vector fields 

d d d 

Y12 = Y 13 = x 2 — , Y M = xv—, 

ox ox ox 

ox Ov ov 

Indeed, using Tables 1 and 2, we readily obtain that [Wq, Vq] C Vq. Moreover, as we already 
know that Wq is a Lie algebra, S(Wg, Vq) becomes a quasi-Lie scheme. 

Since we impose that g+X, which is given by (4.2), must be of the form (8.1) up to a t- 
reparametrisation, we obtain 

11 

g+X = ba(t)Y a = £(t) [(Yi - 3Y 3 - Y 6 ) + f(t)Y 9 + g(t)Y 8 + h(t)(Y 7 + Y 4 )} , (8.2) 

a=l 

for a certain t-dependent non- vanishing function £ (t) . Taking into account that the vector fields 
Yi, . . . , Y11 are linearly independent over K, we see that 

6 2 = 0, 65 = 0, 6io = 0, b n = 0. 

Bearing in mind the form of the coefficients b a given by (4.3) and recalling that a > 0, we see 
that \>2 = entails n = 1. Meanwhile, conditions 65 = and b\\ = entail a = and 

l ir-l {8 - s> 

a at 

respectively. In turn, a = also entails that 69 = 610 = 0. Moreover, from (8.2) we also see that 

63 = — 3&i, b e = -h, 67 = 64. 

From these conditions and n = 1, cr = 0, we obtain 

ao a a^5 a 1 d5 27 1 dao ,^ ^ 

a 5' a 3 5 J 5 dt 5 clq dt 

Since 5, a > 0, it can readily be seen that the second equation is an immediate consequence 
of the first one, which in turn implies ao = —a 2 / 5. Using that 5(0) = a(0) = 1, we obtain 
ao(0) = —1. Previous results along with (8.3) ensure that the last condition in (8.4) holds. 
Hence, 61, 64 and 6s become 

- ao 7 1 dao 3 da 

61 = , 6 4 = ai Tr + --7r> 

a ao at a at 
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a 2 a a\da 2(da\ 1 da 



bs = -— + --77 + hr -—^2, ( 8 - 5 ) 

ao ao at a o a \dt J oq at z 

and the remaining coefficients b a are either zero or can be obtained from them. 

Proposition 8.1. Every Gambier vector field X is a quasi-Lie system relative to the quasi-Lie 
scheme S(Wg, Vq) and the Lie algebra Vq if and only if n = 1, ao(0) = —1 and a = 0. In such 
a case, every t-dependent transformation g € Q(Wg) given by 

a ( da \ , n „. 

x = ax, v = I —x + av I , 



ao V dt 

transforms X into a Lie system 
dx 

— = Sxv -x 3 + b 4 b7 l x + bsbT 1 (x 2 + v) , (8.7) 
ar 

where dr = —aodt/a and b\, 64 and bs are given by (8.5). 



9 Exact solutions for several second-order Gambier equations 

As a result of Proposition 8.1, we can apply to (8.7) the techniques of the theory of Lie systems so 
as to devise, by inverting the change of variables (8.6), new properties of the Gambier equations 
related to second-order Riccati equations. For instance, we obtain new exact solutions of some 
of these families of second-order Gambier equations. 

In the study of every Lie system, special relevance takes the algebraic structure of its Vessiot- 
Guldberg Lie algebras. When a Lie system possesses a solvable Vessiot-Guldberg Lie algebra, 
we can apply several methods to explicitly integrate it (cf. [9, 13, 19]). Otherwise, the general 
solution of the Lie system usually needs to be obtained through approximative methods [51] or 
expressed in terms of solutions of other Lie systems [12]. 

System (8.7) is a special case of a Lie system related to a Vessiot-Guldberg Lie algebra Vq iso- 
morphic to s[(3,M), which is simple and therefore difficult to integrate explicitly. Nevertheless, 
we can prove that (8.7) is in general related to a Vessiot-Guldberg Lie algebra that is a proper 
Lie subalgebra of Vo- Moreover, we next prove that X is related to a solvable Vessiot-Guldberg 
Lie algebra in some particular cases that can be explicitly integrated. 

System (8.7) describes the integral curves of the i-dependent vector field 

X t = Zi + hb^Z 2 + b 8 b^Z 3 . 

For generic t-dependent functions b\, 64 and bs, the above system admits a Vessiot-Guldberg 
Lie algebra spanned by Z\ = X\, Z2 = (A3 + Xj)/2, Z3 = (X% — 2X4) /4 and their successive 
Lie brackets, which generates a proper Lie subalgebra V of Vq. More specifically, we first have 

[Z h Z 2 ] = Z 3 - Z 4 , [Z h Z 3 ] = -{Zi + Z 5 )/2, [Z h Z 5 ] = Z 6 , 

where Z4 = A4, Z5 = A5, Zq = Xq. This shows that the smallest Lie algebra containing Z\, Z 2 
and Z3 must include these vector fields and Z4, Z5 and Z§. Since these vector fields additionally 
satisfy the following commutation relations 

[Z 1 ,Z 4 ] = Z 5 , [Z 3 ,Z 4 ] = 0, [Z 1 ,Z 6 ] = 0, [Z 2 ,Z 3 ] = Z 2 , [Z 2 ,Z 4 ] = -Z 2 , 
[Z 2 , Z5] = Z4 — Z3, [Z 2 , Zb] = Z5 — Z\, [Z 3 , Z 5 ] = (Z5 + Zi)/2, [Z 3 ,Z 6 ] = Z 6 , 
[Z 4 , Z 5 ] = -Z 1} [Z 4 , Z 6 ] = 0, [Z 5 , Z 6 ] = 0, 
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we see that they span a six-dimensional proper Lie subalgebra V of Vq. Moreover, it is easy to 
show that V ~ (Zx + Z 5 , Z 3 - Z 4 , Z 2 ) ©s (Z e , Z x - Z 5 , Z 3 + Z 4 ), where {Z 1 + Z 5 , Z 3 - Z 4 , Z 2 ) ~ 
sl(2,R), namely 

\Z X + Z 5 , Z 3 - Z 4 \ = -2(Zi + Z 5 ), 

[Z 2 , Z 3 - Z 4 ] = 2Z 2 , [Z x + Z 5 , Z 2 ] = 2(Z 3 - Z 4 ), 

the linear space (Zq, Z\ — Z$, Z 3 + Z 4 ) is a solvable ideal of V and ©5 stands for the semidirect 
sum of (Zi + Z 5 , Z 3 - Z 4 , Z 2 ) by (Z 6 , Z\ - Z 5 , Z 3 + Z 4 ). 

Hence, X is related to a non-solvable Vessiot-Guldberg Lie algebra and it is not known 
a general method to explicitly integrate X for arbitrary 61, b 4 and 6s. Nevertheless, we can 
focus on a particular instance of these functions so that X can be related to a solvable Vessiot- 
Guldberg Lie algebra V±. For example, consider the case b 4 = 0. We then have 

Vy = (Z^Z^Zt + Z^Ze). 

Note that the derived series of V\ become zero, i.e. V 1 = [Vi, V\] = {Z\ + Z§,Zq), V 2 = 0. In 
other words, V\ is solvable and we can expect to integrate it through some method. Let us do 
so through the so-called mixed superposition rules, i.e. a generalisation of superposition rules 
describing the general solution of a Lie system in terms of several particular solutions of other 
systems and a set of constants [29]. In our case, it is known (cf. [29, p. 194]) that the general 
solution of (8.7) can be written in the form 

_ MVyiij) + A 2 % j (r) + X 3 v y3 (T) d ( \\v yx {r) + A 2 % 2 (r) + A 3 ^ 3 (r) 

X \ T ) ~ \ „. /_\ 1 \ „ /_\ 1 \ „ /_\ ' V \ T ) 



Aiyi(r) + X 2 y 2 (r) + X 3 y 3 (r) ' dr \ Xiyi(t) + X 2 y 2 {r) + X 3 y 3 (r) 

where (Ai,A2,As) G M 3 — {(0,0,0)} and (yi{T),v yi (T),a yi (T)), with i = 1,2,3, are linearly 
independent solutions of the linear Lie system 

dy dv y da y - - 1 

~=Vy, —=ay, —= bsb, ay. 

The above Lie system is related to a Vessiot-Guldberg Lie algebra spanned by 

d d d d d 

W\ = v y — + a y — , W 2 = a y ^—, W 3 = 2a y — , W 4 = -2a y — , 
ay ov y oa y o v y vy 

which close the same commutation relations as Z\, Z 3 , Z\ + Z5 and Zq. Hence, this system is 
related to a solvable Vessiot-Guldberg Lie algebra and it can easily be integrated: 

a y = exp^J b^{r')b'[ L (T')dT'^j , v y = J a y {T')dr' ', y = j Vy(r')dr' '. 

Since we can assume 

Aiyi(r) + X 2 y 2 {r) + A 3 y 3 (r) = / / exp / ^(r"')^ 1 ^"')^'" \ dr" dr' , 



we obtain the solution of (8.7) for b 4 = 0, i.e. 

r^y{r'h(r'')b- 1 \r'')dr'')dr' = d 

f f' exp (j"" b 8 (T"')bYHT'")dT"^ dT»dT> ' VT dT T ' 



From above results, we have the following proposition. 
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Proposition 9.1. Given a second-order Gambier equation with n = 1, ao(0) = —1, a = and 
a particular solution a(t) with a(0) = 1 of 

1 dan 3 da 

«i = --^--^r, 9.1 
ao at a at 

its general solution reads 

1 J r(<) exp (/ r ' 6 8 (t")&T V)<fr") «V 



a(r(t)) jr(t) jt' exp /jr» 5 8 ( r ///) 6 -i( T ///) dr //A dr // dr / 



s(r(t)) 



where dr = —andt/a and b% and b$ are given by (8.5). 

Note 9.1. Since equation (9.1) can be easily integrated, the above proposition shows that we 
can integrate exactly every second-order Gambier equation with n = 1, ao(0) = —1 and a = 0. 



10 Conclusions 

Two quasi-Lie schemes have been introduced to analyse the second-order Gambier equations. 
Our first quasi-Lie scheme has been used to recover previous results concerning the reduction of 
such equations to reduced canonical forms from a geometric clarifying approach, which allowed 
us to fill a gap in the previous literature. This quasi-Lie scheme also led to determine some quasi- 
Lie systems related to certain second-order Gambier equations, which enabled us to transform 
them into second-order Kummer-Schwarz equations. We have expressed the general solutions of 
such second-order Gambier equations in terms of particular solutions of t-dependent frequency 
harmonic oscillators and Riccati equations. Additionally, new constants of motion were derived 
for some of them. 

The introduction of a second quasi-Lie scheme resulted in the description of an integrable 
family of second-order Gambier equations related to second-order Riccati equations. 
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